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On the basis of a fit to the Sivers effect in deep-inelastic scattering, we make predictions for single- 
spin asymmetries in the Drell- Yan process at BNL RHIC. 
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PACS numbers: 13.88.-l-e, 13.85.Ni, 13.60.-r, 13.85.Qk 

I. INTRODUCTION 

The Sivers effect [1] was originally suggested to ex- 
plain the large single spin asymmetries (SSA) observed 
in p^p ttX (and p^p ttX) at FNAL [2] and recently 
at higher energies in the RHIC experiment [3] . The effect 
considers a non-trivial correlation between (the trans- 
verse component of) the nucleon spin St and intrinsic 
transverse parton momenta pt in the nucleon. It is pro- 
portional to the "(naively or artificially) T-odd" struc- 
ture (St X pT)PAr and quantified in terms of the Sivers 
fimction /f^(a;,p^) [4] whose precise definition in QCD 
was worked out only recently [5-7] . 

A particularly interesting feature of the Sivers function 
(and other "T-odd" distributions) concerns the universal- 
ity property. On the basis of time-reversal arguments it 
was predicted [6] that in semi-inclusive deeply inelas- 
tic scattering (SIDIS) and in the Drell- Yan process (DY) 
have opposite sign (our definition of the Sivers function 
follows the Trento Conventions [8]), 



/it Pt) SIDIS — /itI-^iPt)! 



(1) 



The experimental check of Eq. (1) would provide a thor- 
ough test of our understanding of the Sivers effect within 
QCD. In particular, the experimental verification of (1) 
is a crucial prerequisite for testing the factorization ap- 
proach to the description of processes containing px- 
dependent correlators [9-11]. 

Recent data on SSA from SIDIS [12-14], and in par- 
ticular those from transversely polarized targets [15-17], 
provide first measurements of the Sivers effect in SIDIS. 
On the basis of this information it was shown that the 
Sivers effect leads to sizeable SSA in p^Tr~ — ^ I'^l'X, 
which could be studied at COMPASS [18], and in p^p 
l^l^X or pp^ l^l^X in the planned PAX experiment 
at GSI [19, 20], making the experimental check of Eq. (1) 
feasible and promising [21]. Both experiments, which 
could be performed in the medium or long term, have 
the advantage of being dominated by annihilation of va- 
lence quarks (from p) and valence antiquarks (from p or 
7r~), which yields sizeable counting rates. Moreover, the 
processes are not very sensitive to the Sivers antiquark 



distributions, which arc not constrained by the present 
SIDIS data, see [21-25]. 

On a shorter term the Sivers effect in DY can be 
studied at RHIC in p^p l^l~X. Similar and other 
spin physics prospects at RHIC are discussed in [28-34]. 
Other earlier predictions for SSAs in DY were given in 
[35, 36]. 

In pp-coUisions inevitably antiquark distributions are 
involved, and the counting rates are smaller. In this work 
we shall demonstrate that the Sivers effect SSA in DY 
can nevertheless be measured at RHIC with an accuracy 
sufficient to unambiguously test Eq. (1). In particular, 
by focusing on certain kinematic regions the effect of the 
unknown Sivers antiquark distribution function can be 
minimized. And, by focusing on the opposite kinematic 
regions one can gain first information on the Sivers an- 
tiquark distribution itself. For our estimates we use the 
Sivers function extracted from HERMES data [16] in [22], 
see also [26, 27]. 

It remains to be noted that the theoretical understand- 
ing of SSA in p^p — > nX , which originally motivated the 
introduction of the Sivers effect, is more involved and less 
lucid compared to SIDIS or DY, as here also other mech- 
anisms such as the Collins effect [37] and/or dynamical 
twist-3 effects [38, 39] could generate SSA. Phenomeno- 
logical studies indicate, however, that in a picture based 
on pT-dependent correlators the data [2, 3] can be ex- 
plained in terms of the Sivers effect alone [40-42] with the 
other effects playing a less important role [43, 44]. For 
recent discussions of hadron-hadron collisions with more 
complicated final states (like, e.g., p^p — > jet jet2-^) we 
refer to Refs. [45]. 



II. SIVERS EFFECT IN SIDIS 

The longitudinal SSA in SIDIS observed first [12-14] 
unfortunately cannot be unambiguously interpreted in 
terms of a unique (Sivers [1], Collins [37] or twist-3 [46]) 
effect [47-51], see Ref. [52] for a recent review. All that 
is clear at the present stage is that these SSAs are dom- 
inated by subleading twist effects [53]. The situation 
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changed, however, with the transverse SSA in SIDIS ob- 
served at HERMES and COMPASS [15-17], where the 
Sivers and Colhns effect can be distinguished due to the 
different azimuthal angle distribution of the produced 
hadrons [4]. 

In particular, in the SSA due to the Sivers effect the 
produced hadrons exhibit an azimuthal distribution oc 
sm{(t)—(j)s) in the plane transverse to the beam axis. Here 
(f) and (j)s are the azimuthal angles of respectively, the 
produced hadron and of the target polarization vector, 
with respect to the lepton scattering plane. 

The Sivers SSA measured at HERMES [16] is defined 
as sum over SIDIS events i as 

^^j^^'^-'^^) (2) 

_ Ei - <Ps,i){NH<i>i, 4>s,i) - N^{ct>i,4>s,^ + tt)} 

N'^^^^(j}i,(j)s,i) are the event counts for the respective tar- 
get polarization (corrected for depolarization effects) . In 
order to describe the Sivers SSA defined in Eq. (2) in 
Ref. [22] two major simplifications are made. First, soft 
factors [9-11] are neglected. Second, for the distribution 
of transverse momenta in Di{z,K^) and ,f^{x,p'^) the 
Gaussian model is assumed. 

The Gaussian model certainly oversimplifies the de- 
scription of "unintegrated" distribution or fragmentation 
functions which are an involved issue in QCD [54]. On 
a longer term, an approach to the pT-dependence of the 
Sivers SSA along the lines of the formalism in Ref. [55] 
would be desirable. The Gaussian model, however, pro- 
vides a good effective description of SIDIS and DY data 
within a certain range of low transverse momenta, and 
is sufficient for the purpose of the present work. The 
free parameters, namely the Gaussian widths, are consis- 
tently constrained in [22] by the HERMES data. 

As regards the flavor dependence of the Sivers func- 
tions, there are no strong constraints from the present 
SIDIS data [16, 17]. In fact, in Ref. [23] where this 
has been attempted all fitted distributions but were 
found consistent with zero. In this situation it is ap- 
pealing to invoke additional theoretical constraints. In 
particular, here we use predictions from the QGD limit 
of a large number of colours N^. 

In this limit the nucleon appears as Nc quarks bound 
by a mean field [56], which exhibits certain spin-flavour 
symmetries [57]. By exploring these symmetry properties 
it was proven in a model independent way that in the 
large- iVc hmit [58] 

frri^iPT) — ^firi^^PT) modulo l/Nc corrections, 

(3) 

for not too small and not too large x satisfying xNc = 
0{N^). Analog relations hold for the Sivers antiquark 
distributions.^ Imposing the large- A^c relation (3) as an 



^ For historical correctness we mention that previously (3) was 
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FIG. 1: The u-quark Sivers function ^fyrgj^is^^^ ^ ^ scale 
of about 2.5 GeV^, as obtained from a fit to the HERMES data 
[16]. Shown are the best fit and its 1-a uncertainty. 

additional constraint, and neglecting effects of antiquarks 
and heavier flavours, it is shown [22] that the HERMES 
data [16] can be described by the following 2-parameter 

ansatz and best fit 

/iTsiDis(^) ~ ~/lTsiDis(^) ~ Ax (l — x)' , (4) 

^ -0.17x°-66^1 - , 

with a per degree of freedom of about 0.3. The fit and 
its 1-a uncertainty due to the statistical error of the data 
[16] are shown in Fig. 1. Several comments are in order 
concerning the ansatz and fit result (4). 

The fit (4) refers to a scale of 2.5 GeV^ roughly set 
by the {Q^) in the HERMES experiment [16]. For the 
extraction we used the parameterizations for fi{x) and 
D^{x) from [62, 63] at the corresponding scale. 

Within the Gaussian model, of course, one docs not 
need to work with the "transverse moment" of the Sivers 
function defined as, and in the Gaussian ansatz given by 

^ ydV^/iV(^,P^) (5) 

Gauss (Pt) Siv f±af\ 



where {Px)siv denotes the Gaussian width of the Sivers 
function. However, doing so one benefits from the fact 
that the fit for the moment f^^^°'ix) (in contrast to 
f^{x,p'^)) is nearly insensitive [22] to the value of 
(Pt) Siv "which is poorly constrained by the data and the 
positivity bound [64]. 

The shape of the Sivers function at large x is not con- 
strained by the data [16]. Our ansatz of (1 — x)^ de- 



discussed in the framework of (simple versions of) chiral models 
[59]. But the way in which (3) was obtained there was shown 

to be incorrect [60]. Recently, in Ref. [61] (a more sophisticated 
version of) a chiral model with vector mesons obeying a hid- 
den local flavour symmetry was discussed, in which the Sivers 
function would obey (3). 
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pendence has some theoretical justification [21], although 
there are also arguments [24] for a (1 — x)^ behavior. 

We remark that the result (4) is in good agreement 
with the fit of Ref. [21] to the preliminary HERMES 
data [15] on the Sivers SSA weighted with a power of 
the transverse momentum Phi_ of the produced hadron. 
With such a weighting the SSA can be interpreted (ne- 
glecting soft factors) unambiguously in terms of fix^''"'ix) 
independently of any model for transverse momentum de- 
pendence [4] . Keeping in mind the preliminary status of 
the data [15], this agreement indicates the consistency of 
the Gaussian ansatz within the accuracy of the data. 

The large- iVc motivated ansatz (4) is confirmed by re- 
sults from the COMPASS experiment [17] where a solid 
polarized ^LiD target [65, 66] was used. Neglecting nu- 
clear binding effects and exploring isospin symmetry, in 
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gously for d, q, etc. Thus, the deuterium target is sen- 
sitive to the flavour combination which is suppressed in 
the large- TVc limit, see (3), and for which our ansatz (4) 
yields zero. This is in agreement with the COMPASS 
data [17] showing a Sivers effect from deuterium target 
compatible with zero within error bars. 

The reason why the large Nc picture of the Sivers func- 
tion works at the present stage of art, is due to the fact 
that the current precision of the first data [16, 17] is 
comparable to the theoretical accuracy of the large-iVc 
relation (3). Thus, 1/A^c corrections (and antiquarks ef- 
fects) cannot be resolved within the error bars of the data 
[16, 17]. In future, when the precision of the data will in- 
crease, it will certainly be necessary to refine the ansatz 
(4) to include 1/Nc corrections and antiquark effects. 

The fit (4) described above is to the final HERMES 
data in Ref. [16]. We observe that the fit is compatible 
within the 1-a range with more precise, but preliminary, 
data that has recently become available [67] . Incorporat- 
ing into our fit these data, which are not yet corrected for 
smearing and acceptance effects, would tend to increase 
somewhat the Sivers function, c.f. [24-27]. 

The sign of /{tsidis(^) (^) agreement with 

the physical picture of the Sivers functions discussed in 
Ref. [68]. 

We remark that, in principle, one could include into 
the fit in addition to SIDIS data also the data on SSA 
in the hadronic processes p^p nX or p^p ttX [2, 3]. 
This possibility was not explored in [22] as these SSA 
could also be due to other mechanisms [37 39]. Further- 
more, as shown in [45] it is not clear that factorization 
holds in these processes in terms of pT-dependent corre- 
lators, for reasons that do not apply to SIDIS and to the 
Drell-Yan process. One also has to keep in mind that it 
is a priori not clear whether the leading twist factoriza- 
tion approach is accurate at the (Q^) = 2.5 GeV^ of the 
HERMES experiment as assumed in [22]. Only careful 
analyses (which will include soft factors) of future data 
from experiments performed at different will reveal 
to what extent this assumption is justified. 



III. SIVERS EFFECT IN SSA IN DY AT RHIC 

The process p^p ^ l~^l~X is characterized by the vari- 
ables s = {pi+p2)'^, the invariant mass of the lepton pair 
Q"^ = {ki + fe)^, and the rapidity 



_ 1 ^^ P2 ■ (fcl + ^2) 

^~ 2 ''pi-(A:i + fc2) 



(6) 



where Pi/2 (and ^1/2) indicate the momenta of the in- 
coming proton (and the outgoing lepton) pair. 

Let us consider the azimuthal SSA defined as a sum 
over the events i according to 



^sin(<^-<^s) 



(7) 



where 1,1 denote the transverse polarizations of the 
proton, the polarized proton moves in the positive z- 
direction, and {(p — (t>s) is the azimuthal angle between 
the virtual photon and the polarization vector. Neglect- 
ing again soft factors and assuming the Gaussian model 
for the distribution of transverse momenta, to leading 
order the SSA is given by 



,sm(0-0s) _ o 



^gIuss Eg eg a^i/iTD (3:2) 
Ege^ari/f(a;i)a;2/f(ar2) 



(8) 

where a = u,u, d, d, etc. and the parton momenta xi,2 
are given by a;i,2 = (Q^/s)^/^ e^^. The dependence on 
the Gaussian model is contained in the factor 



"Gauss 



M 



N 



Siv + (PT)unp 



(9) 



If one introduced in the numerator of (7) an additional 

weight qx/Mpf, where denotes the modulus of the 
transverse momentum of the lepton pair with respect 
to the collision axis, the resulting SSA would be inde- 
pendent of a particular model for transverse momentum 
distribution and given by the expression on the right- 
hand-side of (8) but without the factor acluss- 
argued that such a "transverse momentum weighted" 
SSA might be protected against Sudakov dilution ef- 
fects [69]. These effects need not be negligible, consider- 
ing the fact that at RHIC we deal with SSA at consid- 
erably higher scales than in the HERMES experiment: 
(4GeV)2 < g2 < (20GeV)2 typically at RHIC [28] vs. 
(g2) = 2.5 GeV^ at HERMES. 

The DY process in pp-collisions is sensitive to 

./iTdy('^i)./i(-^2) and to /^^ o'yC'^i)/! (■'^2) on an equal 
footing. Thus, even though the effects of the Sivers anti- 
quark distributions seem small and presently not observ- 
able in SIDIS, one cannot expect this to be the case in 
DY at RHIC. 
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FIG. 2: The azimuthal SSA A^'^'^ '^^'^ in DrcU-Yau Icpton pair production, p — > l^l X, as function of y for the kinematics of tlic 
RHIC experiment with -/s = 200 GeV. The left plots show = {4GcY)^, and the right plots = (20GcV)2. The upper plots show 
the effects of varying the Sivers antiquark distributions within the range of model I — see Eqs. (10, 11). The lower plots show instead 
the results of model II. In all cases, the central estimates are based on our fit (4) for the Sivers quark distribution functions from the 
HERMES data [16]. The inner error band (solid lines) shows the 1-a uncertainty of the fit. The outer error bands (dashed lines) show 
the error from varying the Sivers antiquark distribution functions within the ranges specified in Eqs. (10, 11 ), model I for the upper plots, 
and model II for the lower plots. The x-region explored in the HERMES kinematics is indicated. For Q = A GeV we show the estimated 
statistical error for STAR and PHENIX. For Q = 20 GeV the statistical error at STAR and PHENIX is comparable to the asymmetry. 
See text and Table I for a discussion of the planned RHIC II. 



In order to have a rough idea which regions of y (for 
a given Q^) could be more and which regions less sen- 
sitive to Sivers antiquark distribution functions, let us 
introduce two models: 



with 



0.25 = const model I 



e(x) = ± 



model II 



(10) 



(11) 



where for ,f^x^j^^{x) wc will use the result from Eq. (4) 
taking into account the change of sign in Eq. (1). 

Model I is that the Sivers antiquark distribution is 
±25% of the corresponding Sivers quark distribution at 
all X. Model II is that the ratio of the Sivers antiquark 



over quark distribution is fixed by the ratio of the unpo- 
larized distribution functions. The particular ansatz for 

model II ensures compatibility with the large- Nc limit, 
and it preserves the following sum rule [70] — see also 
[21]- 



E 

a=g,Uyd, . 



. 



(12) 



(Note that in the large- A^c limit the gluon Sivers distri- 
bution is suppressed with respect to the quark one [21].) 

Eqs. (10, 11) represent rough models for fi'^^^^{x) 
which are, however, consistent with the HERMES data 
[16], and compatible with the presently known theoretical 
constraints, namely the large- A^c relations (3), positivity 
constraints [64] and the Burkardt-sum rule (12). This 
makes them sufficient for our purposes. 

In order to present estimates for RHIC we strictly 
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speaking should use the fitted Sivers function (and /f (x)) 
evolved to the relevant scale. Instead, let us assume that 



E 2 f-L(l)a fa 



Ep2 fa fa 
a ^aJl Jl 



2 fUl)o- fa 
l^a ^aJlTovJl 



HERMES 



E„2 fa fa 
a ^aJl Jl 



RHIC 



(13) 

It is difficult to quantify exactly the error introduced in 
this way. However, we believe it to be smaller than other 
uncertainties in our analysis. 

When dealing with unintegrated distribution functions 
at higher energies one must take into account a broad- 
ening of the average transverse momentum. Considering 
that this effect due to radiating off soft gluons does not 
depend on polarization, we roughly estimate 



/ 2NRHIC 

xPt/sw 

12 \RHIC 
\PT/unp 



2(p| 



^HERMES 

Siv 

^HERMES 
unp 



0.66 GeV^. (14) 



Similar values for {px)unp at high energies were obtained 
in [42]. Next, let us assume that the transverse moment 
/itdy little affected by that. The combined effect 
of our assumptions is that the SSA decreases at higher 
energies - in qualitative agreement with [69] , though not 
as fast as discussed there. 

We recaU, that from the HERMES data [16] and the 
positivity bound [64] the Gaussian width of the Sivers 
function is only poorly constrained. Under the assump- 
tion (14) one may expect it at RHIC to be in the range 
(Pt)I™^ « (0.20 ... 0.64) GeV^ [22]. However, the only 
place where this value numerically matters is the Gauss 
factor in Eq. (9). Varying (p7-)si™^ the above range 
yields ag^^gg = 0.81 • (1 ± 10%), i.e. it alters our estimates 
for RHIC only within ±10%. 

On the basis of the above assumptions and taking into 
account the change of sign for the Sivers distribution 
function in Eq. (1) we obtain the results shown in Fig. 2. 
We observe that the Drell-Yan SSA is noticeably sensi- 
tive to the Sivers g-distribution functions, as modeled in 
Eq. (10). The effect is more pronounced at larger dilep- 
ton masses Q. However, we note that at Q = 4GeV in 
the region of positive rapidities y the effect is moderate, 
and does not alter the sign of the asymmetry. 

The region of negative rapidities y is strongly sensi- 
tive to the Sivers ^-distribution. For positive e{x) in 
Eqs. (10, 11) the SSA is positive, for negative e{x) it 
is negative. As already mentioned, the effect of Sivers-g 
is more pronounced at larger Q^. Thus, by focusing on 
these regions of y, one could gain the first information 
on the magnitude and sign of the Sivers-g distributions. 

In order to see explicitly in which kinematical region 
our predictions are constrained by data, we show in Fig. 2 
the x-range covered in the HERMES experiment which 
constrained the fit of the Sivers function. In this con- 
text it is worthwhile remarking that our estimates for 
RHIC based on SIDIS data are complementary to those 
made in Ref. [33]. There information on the Sivers func- 
tion was used from the data on SSA in p^p ttX [2]. 
Assuming factorization for this process and considering 



that other mechanisms cannot explain the observed SSA 
at large xp [43, 44], one finds that the data [2] constrain 
the Sivers function at large x > 0.4 above the region 
explored at HERMES. Thus, in the region of large posi- 
tive y, where our estimates are not constrained by SIDIS 
data, see Fig. 2, the estimates of Ref. [33] could be more 
reliable. 



5A 



JLdt 



STAR 



y 4GeV 20 GeV 

-0.5 0.007 0.09 

0.5 0.006 0.06 

1.5 0.007 0.11 



125 pb-^ 



PHENIX 



y 4 GeV 20 GeV 

-1.8 0.008 0.2 

0.0 0.017 0.13 

1.8 0.008 0.2 



125 pb" 



RHIC II 



y 4 GeV 20 GeV 

±2.5 0.003 0.03 

±1.5 0.001 0.01 

±0.5 0.001 0.01 



10 X 125 pb-^ 



TABLE I: StatisticaJ errors SA for the Sivers SSA in Drell Yan for 
the PHENIX and STAR detectors at RHIC: Errors are shown for 
dilepton masses of Q = 4 GeV and 20 GeV assuming an integrated 
luminosity of f Ldt = 125p6~^ and a beam polarization of P = 0.7. 
Error estimates have been carried out using the event generator 
PYTHIA. Projected errors are also shown for a possible future 
dedicated experiment for transverse spin with large acceptance at 
RHIC II (luminosity upgrade); see text for details. 



In Table I we show the statistical errors dA for the 
Sivers SSA in DY estimated with PYTHIA considering 
the acceptance of the STAR and PHENIX detectors. De- 
tector acceptance is conveniently specified in terms of 
pseudo-rapidity 77 = ln(tan|), which is directly related 
to the scattering angle 9 of the lepton pair with respect to 
the beam-axis, and thus to the geometry of the detector. 
In the following acceptance cuts imposed on leptons will 
be given in pseudo rapditity. However, asymmetries and 
their errors are analyzed in bins of photon rapidity y. We 
assume an integrated luminosity J Ldt = 125 pb""'^ and a 
beam polarization of P = 0.7. We use these parame- 
ters as an upper estimate for the statistics these exper- 
iments could acquire with transverse beam polarization 
before RHIC detector and luminosity upgrades will be- 
come available in 2012 (RHIC II). The statistical errors 
can be easily scaled to different parameters for integrated 
luminosity and beam polarization. 

For STAR, which covers the range — 1 < 77 < 2 for 
e+e~, estimates for 6 A are presented for bins centered 
at y = -0.5, 0.5 and 1.5. For PHENIX, which covers 
\ri\ < 0.35 (for e+er) and 1.2 < \ri\ < 2.4 (for ^i+M") we 
have chosen bins centered at y = and \y\ = 1.8. The 
bins in the dilepton mass Q are chosen respectively at 
4 GeV and 20 GeV. 

For Q = 4 GeV the Sivers SSA can be measured at 
STAR and PHENIX. For illustrative purposes the esti- 
mated SA for STAR and PHENIX are shown in Fig. 2. 
The 6A a.t Q = 20 GeV is of the order of magnitude of 
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the asymmetry itself. 

The region of higher Q could, however, be addressed 
taking advantage of the higher luminosity available at 
RHIC II. Wc consider a new large acceptance experiment 
dedicated to Drcll Yan physics with transverse spin. The 
new detector could be loc;atcd at one of the RHIC inter- 
action regions which are not equipped with spin rotator 
magnets and therefore always provide proton-proton col- 
lisions with transverse beam polarization. 

We assume that at RHIC II the luminosities will be 
higher by a factor 2.5 through electron cooling and that 
an additional factor 4 in luminosity will result for the new 
experiment from special focusing magnets close to the 
interaction region. We expect that for a multi-year run 
with the new RHIC II detector an integrated luminosity 
of 10 X 125 pb~^ becomes easily accessible. The proposed 
detector covers an acceptance of \ri\ < 3.0. In Table I we 
present estimates for SA for bins centered at \y\ = 2.5, 
1.5, 0.5. Clearly, RHIC II could access the Sivers effect 
also at large dilepton masses, where the effects of Sivers 
antiquarks are expected to be more pronounced. 

A "symmetric" (with respect to rj) detector will allow 
to gain a factor 2 in statistics for SSA compared to dou- 
ble spin asymmetries, as either beam polarization can be 
used. This additional factor is not reflected in the esti- 
mates presented in Table II. 



IV. gr-DEPENDENCE OF SSA 

In the previous Sections we have discussed SSA inte- 
grated over transverse momenta, namely over qx of the 
lepton pair in the DY process (or, Ph± of the produced 
hadron in SIDIS). However, SSA as functions of qt are 
equally interesting observables. 

In the case of the DY Sivers SSA (8) we were able 
to minimize the model dependence to a certain extent 
(within the Gaussian model), e.g., by showing that the 
estimated SSA varied only within ±10% when the Gaus- 
sian width of the Sivers function was allowed to vary over 
a wide range {pt)s™'^ ~ (0.20 . . . 0.64) GeV^. When dis- 
cussing the DY Sivers SSA as a function of qr our results 
are model-dependent. Nevertheless certain features are 
of a general character and it is worthwhile addressing 
them here. In what follows we assume for illustrative 
purposes that (pI)™^ ~ 0.3 GcV^. 

It is worthwhile stressing that unpolarized DY data on 
the qx dependence can be well described by means of 
the Gaussian model up to qr < (2-3) GeV (depending to 
some extent on the invariant mass Q and the kinemat- 
ics of the respective experiment, sec the detailed study 
in Ref. [42]). (For a description of the qr dependence 
of unpolarized DY in terms of the QCD-based formalism 
of [55] we refer to the work [71] and references therein.) 
It is by no means clear whether the transverse parton 
momentum dependence of the Sivers function can be de- 
scribed equally satisfactory by the Gauss ansatz. This, 
in fact, is among what we are going to learn from RHIC 



and other experiments. 

Let us first mention an unrealistic property of the 
Gauss model. Eq. (8) was obtained upon integrating the 
relevant transverse momentum dependent cross sections 
over qT from to oo. In QCD the corresponding dia- 
grams diverge, and physically it is clear that the large qr 
must be cutoff at a scale ~ 0{Q) [54]. 

Experimentally no artificial large-q'T cutoff is needed, 
since correct kinematics imposes a kinematic limit. What 
is more interesting in our context is to impose a low-^r 
cutoff. This makes sense, and could be even preferable, 
because it could allow to increase the asymmetry, for the 
Sivers SSA has a kinematical zero, i.e. it is SSAoc qr- 

Let us define the ^T-dependent DY Sivers SSA as 



A] 



sin{(t>-(t>s) 



{qr) = 



num(gT) 



'''' ' - den(^Ti 
with the numerator and denominator given by 

^^^^^^^ ((p|)|HIC + (p2^)KHIC)2 
( q'r 



(15) 



unp 



where 

,Siv 



af/i = // dg2dy^^e^^iA-^«;(:ri)x2/f(x2) 



cuts 

auu = jj dQMy^Y.^lx,f^{x,)x2m^2) ■ (17) 

cuts " 



Integrating den (or num) over from to oo yields the 
total unpolarized DY cross section auu (or flQauss'^^B'T')- 

Considering a polarization- independent if factor ~ 1-5 
and under the above discussed assumptions we obtain for 
the (quark dominated) kinematics in the forward region 
1.2 < y < 2.4 for Q = (4-5) GeV the results shown in 
Fig. 3a. (The uncertainty of imm(qT) due to the l-cr 
region of the HERMES fit and the effect of antiquarks is: 
Kf^ctorCTuT = (6.0± 3.5) -10-12 barn. For better visibility 
Fig. 3a shows only the central value.) 

Clearly, by including the low-gr region the SSA is di- 
minished. One could increase the SSA by applying a low- 
qx cut. The effect of such a cut is illustrated in Fig. 3b. 
E.g., choosing g^in = 1 GeV^ increases the SSA by about 
20% which is close to the optimum for our choice of pa- 
rameters. ^ Of course, there is a price to pay for. The 



^ Had we chosen {p^)^™'^ = 0.6 GeV^ the gain would be much 

more spectacular, namely up to a factor of 2.5. The optimum 
would then, however, appear for Qmin = 3 GeV^, i.e. beyond the 
presumed range of applicability of the Gauss model. Also the 
loss of statistics would then be considerable. 
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FIG. 3: a. The numerator and denominator, see Eqs. (15-17), for the azimuthal SSA A^^*^ in the DY process as functions of the 
dilepton transverse momentum qx for 1.2 <y < 2.4 at RHIC with ^ = 200 GeV and = (4GcV)^. The absolute numbers for the cross 

sections are somewliat altered by the assumption (13). b. The azimuthal SSA A^^'^ "^^^ and its estimated statistical uncertainty as 
functions of the Iow-qt cut qmin, J-e. the ratio of the numerator and denominator in Fig. 3a integrated respectively over qx 6 [qmim oo]. 
Both the SSA and its uncertainty are normalized with respect to their values for qmin = 0- 



statistical error grows by ~ 50% 

due to loss of statistics, i.e. fewer counts A''. In princi- 
ple, one can find a favoured qmim which would allow to 
achieve the maximal relative accuracy in the experiment. 
This is a good illustration of the old "golden rule" of 
spin physics: the best analyzing power is in the kinemati- 
cal region where (spin asymmetry)^ x (number of counts) 
reaches a maximal value. 

The above numbers are strongly model dependent. 
The quantitative results change already by choosing dif- 
ferent parameters within the Gaussian model, cf. foot- 
note 2. However, as an important qualitative c;onclusion 
of this exercise, we point out that in principle the intro- 
duction of an appropriately optimized law-qr cut may 
be a helpful device to improve the experimental signal 
for the Sivers SSA in DY. 

Since in this context the crucial feature is the presence 
of a kinematical zero, the same procedure may turn out 
helpful in studies of other SSA or other phenomena re- 
lated to parton transverse momenta in DY, SIDIS and 
other processes. 

V. SIVERS SSA IN DY AT PAX AND COMPASS 

Previously transverse momentum weighted Sivers SSA 
were estimated for the PAX and the COMPASS (hadron 
mode) experiments [21] on the basis of a Sivers func- 
tion fit to the preliminary HERMES data on the Ph±- 
weighted Sivers SSA [15]. Here, we solidify the estimates 
and conclusions of Ref. [21] on the basis of the Sivers 
function fit (4) to the Gnal HERMES data [16]. 

The fixed target mode at PAX would make available 
s = 45GeV^, which would allow to access lepton pairs 
in p^p l~^l~X in the region = 2.5 GeV^ below 
the J/-!/)-resonance and well above the region of dileptons 



from ^(1020) decays. This corresponds to the average 
scale of the HERMES experiment [16], thus at PAX it is 
unnecessary to consider pr-broadening effects as we have 
modeled in Eq. (14). 

On the basis of the fit (4) we obtain for the PAX exper- 
iment the estimates plotted in Fig. 4. Nearly the entire 
range of y for this kinematics is constrained by results 
from HERMES data. The SSA is sizeable, (5-10)% in the 
central region, and could be measured, see [19]. The un- 
certainty due to the unknown Sivers antiquark distribu- 
tion functions is completely negligible, as expected [21]. 

In the COMPASS experiment, using the possibility of 
a hadron beam, one could study SSA in p^7r~ — > X 
at s = 400 GeV^ and, e.g., ^ 20GeV^ well above the 
resonance region of J /ip and other charmonia. In this 
kinematics pr-broadening effects can be expected to be 
comparable to RHIC, and we use estimates analogous to 
(14). 

Our estimate for the COMPASS experiment is shown 
in Fig. 4. Again we observe that for this kinematics one 

probes the Sivers function in the x-region, where it is well 
constrained by the HERMES data [16]. Also at COM- 
PASS, the Sivers SSA is sizeable, about (4 - 8)% in the 
central region, i.e. sufficiently large to allow a test of the 
prediction (1). In order to see this, let us make the fol- 
lowing crude estimate. The differential cross section for 
-K~p ji'^ ix~ X is about da/dQ ~ 50pb/GeV in the 
region y > in a kinematics comparable to that con- 
sidered here [72]. Assuming an integrated luminosity of 
lO^^cm"^, which is realistic [18], one could measure with 
an ideal detector and a target polarization of P = 100% 
the SSA with an accuracy of 5 A = l/(P\/iV) 0.5% in 
the forward region y > in a bin of size AQ ~ 1 GeV. Of 
course, this accuracy cannot be achieved with a realistic 
detector and polarization — but this is not necessary for 
a first test of the change of sign in Eq. (1). 
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FIG. 4: The azimuthal SSA "^^^ as function of y in Drcll-Yan Icpton pair production in p^p i+Z X for the kinematics of the 

PAX fixed target experiment with s = 45 GeV^ and = 2.5 GeV'^, and in Drell-Yan lepton pair production in p^n^ —> I'^l^X for the 
kinematics of the COMPASS experiment with s = 400 GeV^ and = 20GeV^, respectively. The estimates are based on the fit for 
the Sivers q-distribution functions, see Eq. (4), obtained from the HERMES data [16]. The inner error band (sohd hnes) shows the 1-a 
uncertainty of the fit. The outer error band (dashed hnes) arises from assuming that the Sivers q-distribution functions are proportional 
to the unpolari2sed antiquarks see Eqs. (10, 11). For the PAX experiment the uncertainty due to Sivers antiquarks is not visible on this 
scale. The x-region explored in the HERMES kinematics is shown. 



The impact of Sivers antiquark distributions is weak 
and negligible as noted in Ref. [21]. In Fig. 4 we demon- 
strate this for the model II in Eq. (11). For the estimates 
in Fig. 4 wc used the paranicterizations [62, 73] (using 
[74, 75] confirms the results). 

Thus, though at the present stage one cannot provide 
dedicated estimates of SA for these experiments as wc 
did for RHIC in Sec. Ill, we conclude that also PAX and 
COMPASS could test the prediction (1) for quarks. How- 
ever, both experiments are insensitive to Sivers antiquark 
distributions [21], which underlines the unique feature of 
RHIC with respect to this point. 

We stress that quantitative analyses and extractions 
of the Sivers function from COMPASS, PAX and RHIC 
data will require a good theoretical control on NLO- 
QCD-corrections, soft factors and evolution effects. For 
discussions of some of these issues in the context of 
(collinear) double spin asymmetries in particular for the 
PAX kinematics see [76, 77]. 

VI. CONCLUSIONS 

We conclude that PHENIX and STAR can confirm the 
change of sign of the Sivers ^'-distribution function in 
SIDIS and the Drell-Yan process in Eq. (1) (as can PAX 
with an antiproton and COMPASS with a pion beam). 



Both PHENIX and STAR can provide first information 
on the Sivers g-distribution functions which are not con- 
strained by the current SIDIS data from HERMES or 
COMPASS, unUke PAX and COMPASS. 

There are more recent (preliminary!) HERMES data 
[67]. Our fit is compatible with the new data, but the 
tendency is to increase the Sivers fimction. This tendency 
would result in even more optimistic estimates for DY at 
RHIC (and COMPASS and PAX). 

We confirm Refs. [24, 25] with respect to the capabil- 
ity of RHIC to access the Sivers g-distribution. In ad- 
dition we point out the possibility to learn about Sivers 
^-distribution from RHIC. 
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